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Abstract 

We review the recent theoretical progress in the construction and solution of the evolution 
equations which govern the scale dependence for the twist-three structure and fragmentation 
functions of the nucleon. 



1 Introduction. 



In the naive parton model the cross sections for several inclusive processes can be expressed 
in terms of the density of probability of partons (0) in a hadron (h) F^/h(x) or hadrons in 
a parton Thu{C) with x, l/( being the momentum fraction of the final state particle with 
respect to the initial one. The probabilistic picture relies on the fact that the constituents 
in the high energy processes behave as a bunch of noninteracting quanta at small space-time 
separations. However, the description of hard reactions using this simple intuitive picture 
is very restrictive. The rigorous field theoretical basis at the leading twist and beyond 
comes from the asymptotically free QCD and the use of the factorization theorems which 
give the possibility to separate the contributions responsible for physics of large and small 
distances involved in any hard reaction. At the lowest twist level, the parton model is 
trustworthy and can be mapped onto the language of operator product expansion approach 
(OPE) (for deep inelastic scattering). The contribution of large distances is parametrized 
by the distribution (fragmentation) functions mentioned above, or by parton correlators in a 
broader sense, which are uncalculable at the moment from the first principles of the theory. 
The second ones — hard-scattering subprocesses — can be dealt perturbatively. The parton 
distributions are defined in QCD by the target matrix elements of the light-cone correlators 
of field operators . This representation allows for the estimation of these quantities by the 
non-perturbative methods presently available: like the QCD sum rules |J, the effective 
chiral Lagrangians §, the MIT bag model ||,§||@§g. 

After addressing for over decades mainly the spin averaged observables which have probed 
the nonpolarized substructure of the hadrons and have yielded an important information on 
the parton content of nucleon, the attention has been shifted towards more subtle dynamics 
underlying the polarized scattering. A renewed interest in the high energy spin physics in 
the last years has been concentrated on the transverse spin phenomena in hard processes. 
It revives many ideas developed over a decade ago. In particular, the notion of the twist-2 
transversity distribution hi (x) , first mentioned by Ralston and Soper |]15| has been reinvented 
as well as its twist-3 counterparts have been addressed fT6|. Due to chirality conservation 
h\(x) cannot appear in the inclusive deep inelastic scattering (DIS) but it can be measured, 
for instance, in Drell-Yan reactions through the collision of the transversely polarized hadrons 
13, TS] and in semi-inclusive pion production in DIS on the nucleon ||17|| . In the last case, 



it enters into the cross section as a leading contribution together with the twist-3 chiral-odd 
spin-independent fragmentation function T{Q fTT|| . It is well known that there is considerable 
difference between the structure and fragmentation functions. Namely, the moments of the 
former are expressed in terms of reduced matrix elements of the tower of the local operators 
of definite twist. This property is established by exploiting the Wilson OPE for inclusive 
DIS. Although the light-cone expansion for the fragmentation processes is similar to DIS, 
the moments of the corresponding functions are not related to any short-distance limit. As 
a substitute for the local operators come the Mueller's time-like cut vertices |18|], which are 
essentially nonlocal in the coordinate space so that the analogy to the operator language is 
only useful mnemonic. 

The transverse spin phenomena |TS[ attained by the inclusive DIS are associated with 



explicitly interaction-dependent, i.e. higher twist (twist-3), effects. To disentangle the un- 
derlying complicated dynamics, the unravelling of the twist-3 effects in the hard processes, 
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which manifest the quantum mechanical interference of partons in the interacting hadrons 
is needed. At the twist-3 level, the nucleon has three structure functions g2(x), e(x) and 
hi{x) and the fragmentation functions corresponding to each introduced distribution. These 
new characteristics have been the subject of intensive theoretical study until recently since 
they open a new window to explore the nucleon content. The most important advantage 
of the twist-3 structure functions is that while being important for understanding of the 
long-range quark-gluon dynamics they contribute at the leading order in 1/Q (Q being the 
momentum of the probe) to certain asymmetries || [IU| and, therefore, are directly accessible 



by either the polarized (semi) inclusive DIS |2(J or Drell-Yan or hadron production in e + e" 
annihilation processes. To confront the theory with high-precision data, the knowledge of 
the size of the logarithmic violation of the Bjorken scaling by the QCD radiative corrections 
in the measurable quantities is highly required. 

There exist two equivalent (and complementary) approaches to studying the (^-depen- 
dence of the structure functions in the leading logarithmic approximation (LLA) of the 
perturbation theory. The first is based on the use of the OPE for the product of currents. 
It makes possible the study of logarithmic violation of Bjorken scaling as well as of the 
power suppressed contributions (higher twists) responsible for many subtle phenomena in a 
polarized scattering. The former is achieved by the calculation of the anomalous dimensions 
of the corresponding local operators. However, there exists an alternative approach to the 
analysis of the corresponding quantities which is based on the evolution equations [pi], Q . In 
spite of the fact that the latter has some difficulties as compared to the former, in the study 
of higher twists (like the loose of the explicit gauge and Lorentz invariance of calculations 
and also the presence of the overcomplete set of correlation functions) it has an important 
advantage as being the closest to the intuitive physical parton-like picture. There is another 
advantage of the latter approach to studying the higher twist effects from the point of view 
of experimental capabilities since the OPE provides us with the moments of the structure 
functions, and in order to extract the former, one needs to measure the latter in the whole 
region of the momentum fraction very accurately. Obviously, it is a quite difficult task even 
for the next generation of colliders. While, with a set of evolution equations at hand, one 
can find, in principle, the Q 2 - dependence of the cross section in question by putting the 
experimental cuts on the region of the attained momentum fractions. This approach can 
also be used in the situations when the OPE is no longer valid, i.e. the time-like processes. 

Thus, the Q 2 -evolution of the parton distributions [2T[ Q can be predicted unambiguous- 
ly by exploiting the powerful methods of renormalization group (RG) and QCD perturbation 
theory. As distinguished from the leading twist evolution, the twist-3 two-quark fragmen- 
tation functions receive contribution from the quark-gluon correlators even in the limit of 
asymptotically large momentum transfer. To solve the problem, one should correctly ac- 
count for the mixing of correlators of the same twist and quantum numbers in the course of 
renormalization. 

In the subsequent discussion we attempt to review the recent theoretical progress in the 
study of the twist-3 polarized and nonpolarized chiral-odd and -even structure and fragmen- 
tation functions in the framework of QCD. The presentation will be organized as follows: 
The first chapter is devoted to the consideration of the parton distribution functions. Here 
we address the issues of the sum rules, construction of the basis of independent correlators 
closed with respect to renormalization group evolution, the machinery for the calculation of 
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the evolution kernels, simplification which occurs in the limit of a large number of colours 
and the orbital momenta. The last features make possible the finding of the analytical so- 
lution of the approximated evolution equations. The second chapter concerns the case of 
fragmentation functions. We mainly cover the same subjects as in the first part of the paper. 

2 Space-like cut vertices. 

2.1 Correlation functions of nonleading twist. 

As we have mentioned in the introduction, the parton distribution functions in QCD are 
defined by the Fourier transforms along the null-plane of the forward matrix element of the 
parton field operators product separated by an interval A on the light cone (n 2 = 0): 

JF(A) = JF(A, 0) = 0*(O)$[O, \n](j)(\n), (1) 

where <fi denotes a quark ip or a gluon field B p and $ is a path ordered exponential along 
the straight line which insures the gauge invariance of the parton distribution 

$[x,y] =Pexp^(x-y) M ^da5 M (y + a(x-y))) . (2) 

We suppress the dependence on the renormalization scale fiR in Eq. (|l|), necessary to render 
this quantity well defined in the field theory. The Fourier transformations from the coordinate 
to the momentum space and vice versa are given by 

F(x) = J ^e lXx (h\F(X)\h), (h\F(\)\h) = J dxe- iXx F(x). (3) 

Both of these representations display the complementary aspects of the factorization. The 
light-cone position representation is suitable to make contact with the operator product 
expansion approach, while the light-cone fraction representation is appropriate for establi- 
shing the language of the parton model. Throughout the paper we will use the light-cone 
position and the light-cone fraction representations in parallel. 

It is well known that in order to endow the theory with parton-like interpretation and to 
get much deeper insight into the corresponding perturbative calculations, it is necessary to 
use to ghost-free gauges. Owing to this fact we choose in what follows the light-cone gauge 
B + = n^B^ = for the boson field. The advantage of this gauge is that the gluon field 
operator B p is related to the field strength tensor G pa via the simple relation 

1 r°° 

B^Xn) = d-'G^Xn) = - / dze{X - z)G +fl (z), . (4) 

Z J —CO 

Here, the residual gauge degrees of freedom are fixed by imposing antisymmetric boundary 
conditions on the field, which allows a unique inversion. Thus, the gauge invariant result 
can be restored after all required calculations have been performed. 

To trace the origin of the operator definition of the hadron's parton density we sketch 
briefly below the factorization procedure of Ellis-Furmanski-Petronzio (EFP) |23| for the 
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hadron matrix element of the T-product of the electromagnetic currents T(P, q) whose imag- 
inary part defines the familiar structure functions of DIS. To this end we use the Sudakov 
decomposition of the four-momentum of the active parton in transverse and longitudinal 
components 

F = xp" + an» + fc£. (5) 

Here n is a light-cone vector n 2 = normalized with respect to the four- vector P = p+ ^M 2 n 
of the parent hadron h of mass M, i.e. nP = 1, and p is a null vector along the opposite 
tangent to the light cone such that p 2 = 0, np = 0. 

The physical observable T(P, q) can be factorized into the hard H and soft S blocks (we 
omit the Lorentz indices of the currents J M ) 

T(P,q)=i J d 4 ze^(h\T{J(z)J(0)}\h) = J \{d%[H{h, q 2 )S(k t ,p, A 2 )] 

= J Y[d%dxiS(xi - (kin))[(H(xip, q 2 ) + . . .)S(k h p, A 2 )] 

i 

= Jl[dx l [(H(x t p,q 2 ) + ...)S(\ u p,A 2 )] = [H0S]. (6) 

i 

where we have used the collinear expansion of the momentum of the struck parton with 
respect to the large +- direction of the hadron momentum. 

S(xi,p, A 2 ) = J n^ 4 M(^ - (h n ))S(ki,p, A 2 ) (7) 

i 

with 

S(k,p,A 2 ) = J nrfV Ei " ifc< (A|0(^l)0(^)...^n)|/l). (8) 

i 

Simple manipulations allow one to write the final answer for the soft part 

S(x u p,A 2 ) = |n^ eEiiAiXi ^l0(Airi)0(A 2 n)...0(A n n)|/i). (9) 

Note that $ = 1 in the gauge we have chosen. By exploiting the Poincare invariance of 
the forward matrix element we can exclude the overall translation and, in a particular case, 
come to Eq. ([J). 

The multiparton distributions corresponding to the interference of higher Fock compo- 
nents in the hadron wave functions that emerge at the twist-3 level are the generalizations 
of (|T]) to the 3-parton fields and present already in (Q) 

JF(A, p) = F{\ 0, fi) = 0*Gun)0(O)0(An). (10) 

We do not display the quantum numbers of the field operators since they are not of relevance 
at the moment. The direct and inverse Fourier transforms are 

F(x,x') = [ ^^e ax -^>|^(A,/i)|/i), (h\F(\, (i)\h) = I dxdx'e- iXx+ifiX 'F(x,x r ). 

J 2,71 Z71 J 

(11) 
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The variables x and x' are the momentum fractions of incoming and outgoing cff partons, 
respectively. The restrictions on their physically allowed values come from the support 



properties of the multiparton distribution functions discussed at length in Ref. |[24|| , namely 
F(x,x') vanishes unless < x < 1, < x' < 1. 

Beyond the leading-twist level the intuitive parton-like picture is not so immediate, as 
one usually starts with an overcomplete set of correlation functions. However, the point is 
that the equations of motion for field operators imply several relations between correlators, 
and the problem of construction of a simpler operator basis is reduced to an appropriate 
exploitation of these equalities. The guiding line to disentangle the twist structure is clearly 
seen in the light-cone formalism of Kogut and Soper P5| |. Consider, for instance, the cor- 



relators containing two quarks ipip. Then, decomposing the Dirac field into "good" and 
"bad" components with Hermitian projection operators V± = |7=p7± : i>± = we have 

three possible combinations ip+ip+, ip+ip- ± ipL.ip+, and tp_ip-, which are of twist 2, 3 and 
4, respectively. The origin of this counting lies in the dynamical dependence of the "bad" 
components of the Dirac fermions 

^- = - i -d- 1 (tfl ± + m) 1+ ^ + . (12) 

These components depend on the underlying QCD dynamics, i.e. they implicitly involve 
extra partons and thus correspond to the generalized off-shell partons, which carry the 
transverse momentum. For this reason we come back to the on-shell massless collinear 
partons of the naive parton model, but supplemented with multiparton correlations through 
the constraint (fl2]) . The operators constructed of the "good" components only were named 
quasi-partonic [26]. The advantage of handling them is that they endow the theory with a 
parton-like interpretation for higher twists. 

The EFP approach is close to the OPE; this equivalence is established by identifying the 
moments of the parton correlation functions (S"-block) with the reduced matrix elements of 
the local operators. The singularities of the product of the currents on the light-cone are 
absorbed in the coefficient functions H in front of operators. In the momentum space they 
result in the inverse powers of the large momentum scale Q at which the operators contribute 
to the cross section and it is controlled by their twist r (r = d c — s, where dc is a canonical 
dimension of an operator and s is its Lorentz spin). The leading contribution comes from 
the operators of twist r = 2: 

[<f (O)0(z)] tw - 2 = ]T ...z, n {<f (0)<9/A • • • 9^(0) - traces} . (13) 



However, as has been established in Refs. [27, one can give the definition of the twist 



without appealing to the concept of the local operators which is particularly useful in cases 
when the short distance expansion is no longer relevant, as it happens, for instance, for 
inclusive production of hadron in the e + e~-annihilation et ctr. The point is that the nonlocal 
string operator given by Eq. (p~3|) obeys the Laplace equation 

□ [0*(O)0(;z)] tw ~ 2 = 0, here □ = d 2 (14) 

with the boundary condition on the light-cone 

<f (O)0(z) = £ 1 ...z, n {<P(0)d m d n . . . 8^(0)} > for z 2 = 0. (15) 
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The solution can be written in the form 



oo i / 2\™ i 

[^(O)0(z)r- 2 = r(O)^) + E r, -pn -T / ^(^WH, 

n=0 i J- \ 4 / J0 

(16) 

here d is a dimension of space-time and v = 1 — v . For the time-like processes the connection 
to the local operators is lost, so that we are left with (|T6| ) up to an arbitrary solution of (|T^) 
which vanishes on the light cone. To this accuracy we can define 

00 1 / 7 2 \ n r°° 

(17) 

where the integration region mimics the physical domain of the parton momentum fraction 
of the annihilation channel and thus the matrix elements of this string operator possesses 
the correct support properties in the light-cone variables. 

2.2 g 2 {x). 

Recently, the first experimental data for the measurement of the transverse spin structure 
function g2{x) in the deep inelastic scattering of the longitudinally polarized muon beam on 
the transversely polarized proton target have been reported |20|| . Although at present the 
statistics is too low to be able to extract its perturbative evolution it proves to be important 
to know the theoretical prediction for the latter from QCD. This issue has been extensively 
studied in the literature, therefore, we just outline the main results referring the interested 



reader to the original works p9, [3C], |31, |27], |32], |33[ and reviews p4 



The function we are interested in appears as coefficient in the Lorentz decomposition of 
the antisymmetric part of the hadronic tensor, relevant to the polarized scattering, over the 
appropriate tensor structures: 



i-Im i J dze i( - qz \h\T{j { ^z)J u] {Q)} \h) 



-^^upaq P { S CT #l(x, Q 2 ) + g 2 (x, Q 2 ) (s a - J^zVa 



Using the nonlocal light-cone OPE we can express, as has been discussed in the preceding 
section, the polarized structure functions in terms of the hadronic matrix elements of Fourier 
transformed string operators 

S+9i{x) = y^e^(h\m7 + lB^n)\h), 

S^9t(x) = \j^e^(h\miH^n)\h), (19) 

where S^r denotes the transverse polarization vector of the hadron h (S 2 = —M 2 ) and 
9t = 9i + 92- 

As we have noted above, the higher-twist two-quark operators mix with multiparton 
correlators. Moreover, the operator corresponding to does not possess a definite twist, 
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and as a consequence could not be renormalized multiplicatively. Taking into account the 
equation of motion for the quark field and equality arising from the use of the Lorentz 
invariance, one can find [001 |34| 



xg T (x) - M{x) - K(x) - J dx' D(x, x') = 0, (20) 

/x f \ °~ - r ,D(x,x')+ D(x,x') ,_,„. 

xgAx) = xgrlx) — x—K(x) — x ax r . (21) 

ox J \x' — x) 

Here, we have introduced the new correlation functions 

S^M(x) = l -J ^e iXx (h\mrnj+lH^n)\h), (22) 
S^K(x) = ~ J ^%h\m 1+ d^{\n)\h), (23) 

S^D 1 (x,x') = \J ^^e^-^>|Vi(HgT + ^(0)7 CT ± T5V'(An)|M, (24) 
S^D 2 {x\x) = I J ^^e^'-^(/i|VS(An)g7 + 7^ ± (0)75^(H|/i), (25) 

D(x,x') = l - [D x {x,x') + D 2 (x',x)] (26) 

is a C-even combination of correlators which can enter into the cross section due to even 
photon state in the t-channel under the charge conjugation . The derivative in the correlation 
function K(x) acts on the quark field before setting its argument on the light cone. 

Solving the system of the differential equations (0) and ( plj) with respect to gr(x) the 
integration constant can be found from the support properties of the distribution: gr{x) = 
for |x| > 1. The solution provides us with the following relation between these functions: 



and 



9t(x) = [j9i{P) + \M{x)- [ d -^M{P) 



+ 

where 



d/3 f d(3' 



(3 J P'-p 



(27) 



Y(x,x') = (x-x')D(x,x). (28) 



Here Y(x, x') is explicitly gauge invariant distribution since D is gauge variant provided we 
use a gauge other than the light-cone. To see this, we exploit the advantages of the light-cone 
gauge, where the gluon field is expressed in terms of the field strength tensor by Eq. (|J) and 
take into account the relation 

1 r d\ ±!Ax _ ± J_ PV I e ±^ (29) 

2 J 2tt v ; 2tt x K J 

we can easily obtain the expressions of the gauge-invariant quantities in terms of three- 
particle string operators which are nonlocal generalization of the Shuryak-Vainstein operators 
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z S a [E9|. Generically 



SiYfaJ) = If^e^-^'ihl+S^O,^), (30) 

Z J 271 

S^Y 2 (x',x) = i f^e^-^ihl-S^Am), (31) 

Z J Z7T 



with 



S„(A,0,//) = ^(/zn)*g 7+ zG^ + (0)±7 5 G*^ + (0) ^(An), (32) 



where = |e Ml//0<7 G pCT is the dual field strength tensor and we have used the relation 
ej; a Gp + = with the two-dimensional antisymmetric tensor e^ CT = e p CT +-- 

Thus, the leading order analysis (|27j ) suggests that the structure function can be written 
as the following sum: 

g 2 (x)=g 2 vw (x)+g 2 (x), (33) 

where 

g 2 vw (x) = -g 1 (x) + [jgm (34) 

is the twist-2 Wandzura-Wilczek contribution to the structure functions |3l|, while g 2 (x) is 
a genuine twist-3 (explicitly interaction-dependent up to unessential quark-mass kinematical 
contribution) part which is expressed via the integral of the matrix element of the nonlocal 
operators which measure the quark-gluon correlation function in the target nucleon. 

The distribution functions defined by Eqs. (]T9|), (|2^)-(p5|) form the redundant basis 



of operators closed under renormalization group evolution. Going over to the operators 
composed of the "good" components only, we are forced to consider the renormalization of 
the correlators[]M(a;) and Y(x, x'). 

2.3 e(x). 

In the unpolarized case we define the following redundant basis of the chiral-odd twist-3 
correlations: 

e{x) = X -J^ x %h\m^n)\h) ) (35) 
M(x) = \f ^e iX *(h\mrni + ^n)\h), (36) 

Z J ZTl 

Di(x,J) = \J ^e^-^\h\^n)g 1+ ^(0)^(Xn)\h), (37) 

D 2 (x',x) = \f ^e^'-^{h\^{\n)g^\()) 1+ ^n)\h). (38) 
Z J Ztt Ztt 

The functions D 1 and D 2 are related by complex conjugation [Di(x,x')]* = D 2 (x',x). The 
quantities determined by these equations form a closed set under renormalization; however, 



1 In this discussion, we restrict ourselves to consideration of the non-singlet channel only. 



s 



they are not independent, since there is a relation between them due to the equation of 
motion for the Heisenberg fermion field operator: 

e(x) - M{x) - J dx'D(x, x') = 0, (39) 
where again we have introduced the convention 

D(x, x') = - [D^x, x') + D 2 (x', x)} . (40) 
This function is real- valued and antisymmetric with respect to the exchange of its arguments: 
[D(x, a/)]* = D(x, x'), D(x, x') = -D(x', x). (41) 

Below, in section |2.8| , as an illustration of the self-consistency of the whole approach to 
the higher twists we present a set of coupled RG equations for the correlation functions 
determined by Eqs. fl3"5])-(|3"ED derived in the abelian gauge theory. Relation ( |3"9"| ) provides a 
strong check of our calculations^. It allows the reduction, as we have mentioned above, of 
the RG analysis to the study of scale dependence of the three-parton D and mass- dependent 
M correlators only. 

Introducing the quantity 

Z(x, x) = (x- x')D(x, x'). (42) 

we can easily obtain from Eqs. tfHty and flBlf ) the definition of the gauge-invariant quantities 
Z in terms of three-particle string operators, namely 

Z(x,x') = l -J ^e^'^' (h\Z(X^) + Z(-im, -X)\h), (43) 

where 

Z(X, /i) = Z(X, 0, /i) = ^n)gG +p (0)a^(Xn). (44) 
In the same way, for a mass- dependent non-local string operator 

777 _ 

M 3 (X) = M 3 (A, 0) = -i;(0)l+(iD + (X)) 3 ^(Xn), (45) 
the Fourier transform is 

r rlX 

M j (x) = x j M(x) = / —e lXx (h\M j {X)\h). (46) 

For the spin-dependent scattering discussed below, the only difference is that one should 
insert also a 75-matrix between the fields in the definitions of the string operators Q4"4]), (|4~5]). 

2 This fact follows from general rcnormalization properties of gauge- invariant operators as one expects that 
the counter term for the equation of motion operator can be given only by the operator itself. Its matrix 
element, being taken with respect to the physical state, decouples completely from the renormalization group 
evolution. 
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2.4 h L (x). 

Analogously, the set of correlation functions for the polarized case is as follows: 



h(x) = \S^J ^e^(h\mi<l^n)\h), (47) 

h L {x) = |/ ^e^(h\m^ + -^(Xn)\h), (48) 

M(x) = l -f ^e iXx (h\m^l+l^(Xn)\h), (49) 

K(x) = \f ^e iXx (h\mil + @±l^(\n)\h), (50) 

Dfaa>) = I / ^% lXx -^\h\^n)g 1+ $\0) l5 4,(\n)\h), (51) 



D 2 (x',x) = lj ^ge^'-^(/ i |^(An)g 7+ ^(0)75^(H|/ i ), (52) 
Besides the identity arising from the equation of motion 

h L (x) - M{x) - K{x) - J dx'D(x, x) = 0, (53) 
there is an equation provided by the Lorentz invariance 

2xh 1 (x) = 2h L (x) - x^-K(x) -2x\ dx' ^p^- . (54) 

ox J (x' — X) 

It means that both parts of this equality are expressed in terms of matrix elements of different 
components of one and the same twist-2 tensor operator, and thus should possess the same 
anomalous dimensions. Again, we have introduced the C-even quantity D, which has the 
properties 

D(x,x') = D(x,x'), D(x, x) = D(x', x). (55) 
Following the same line as above, we come to the equation 



d P, to s r?/ n o 2 f ld P 



r 1 d(3 


f dp' 


' d 


d ' 


L (3 2 J 






dj'_ 



h L (x) = 2x 2 I -^h 1 ((3) + M(x)-2x 2 ^ ^M(P) 

Z(f3,f3'). (56) 

A similar relation was found by Jaffe and Ji f\ in Ref. JTOf. Here the dynamical twist- 
3 contribution is explicitly related to a particular integral of the three-parton correlation 
function Z. In terms of local operators it looks like 

n 

(n + 3)[h L ] n = 2[/i 1 ] n+1 + (n + l)M n + - I + l)Z l n , (57) 

i=i 

and the definition of moments of distribution functions is given by Eq. (|112|) . 

As before, excluding the functions fl48|) and (p0|) , and using the relations fl53|) and 
we can chose the basis of independent functions in the form: hi(x), M(x), D(x,x'). 

3 The corresponding expressions in Ref. [JlO] contain misprints. 



10 



2.5 Construction of the evolution equations. 



As long as z 2 7^ 0, the renormalization of the T-product of the operators T {<p*(0)(f)(z)} 
is trivial and reduced to the familiar renormalization constants Z of the fundamental field 
operators entering into the Lagrangian density. However, if z = or z 2 = an additional 
divergence enters the game. This is a well-known fact from the renormalization theory since 
the product of (at least) two field operators entering into the same space-time point (or 
on the light cone) produces an ill-defined quantity from the point of view of the theory of 
distributions and the corresponding infinities have to be regularized and subtracted. In the 
momentum space this results in ultraviolet (UV) divergences of the momentum integrals in 
the perturbation theory, and as a by-product this causes the logarithmic dependence of the 
parton densities on the normalization point. This dependence is governed by the renormaliza- 
tion group. The evolution equations for the leading twist correlation functions determining 
their Q 2 dependence can be interpreted in terms of the kinetic equilibrium of partons inside 
a hadron (for distribution functions) or hadrons inside a parton (for fragmentation function) 
under the variation of the ultraviolet transverse momentum cut-off ||21]| . However, beyond 
the leading twist the probabilistic picture is lost due to a quantum mechanical interference 
and more general quantities emerge, i.e. multiparticle parton correlation functions, whose 
scale dependence is determined by the Faddeev-type evolution equation with a pair-wise 



particle interaction 26 . 
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Figure 1: Ladder diagram for deep inelastic scattering. 

There are two sources of the logarithmic dependence of the correlation functions. The 
first is the divergences of the transverse momentum integrals of the particles interacting with 
the vertex and forming the perturbative loop. Another source is the divergences due to the 
virtual radiative corrections. In the renormalizable field theory the latter are factorized into 
the renormalization constants Z of the corresponding Green functions. However, owing to 
the specific features of the renormalization in the light-like gauge, extensively reviewed in 
the next section, there is mixing of correlation functions due to the renormalization of field 
operators. The latter fact is closely related to the matrix nature of renormalization constants 
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of the elementary Green functions in the axial gauge. For example, after renormalization 
of the fermionic propagator the matrix structure of the bare vertex could be changed, in 
general, since the renormalization matrix acts on the spinor indices of the vertex (see Eq. 

dH)). 

In the Leading Logarithmic Approximation (LLA) there is a strong ordering [|7j of 
transverse particle momenta as well as their minus components, so that only the particles 
entering into the divergent virtual block S, T, II or the particles adjacent to the vertex can 
achieve the maximum values of \k±\ and a (see Fig. [I]): 

\k n ±\ < ••• < IM < |fcuj < A, 

a n ... a 2 <^ ocx, (58) 

while the plus components of the parton momenta are of the same order of magnitude 

x n ~ ... ~ x 2 ~ X\ (59) 

for a n-rank ladder-type diagram. 



2.6 Renormalization in the light-cone gauge. 

A peculiar feature of the light-like gauge is the presence of the spurious IR pole l/k + in the 
density matrix of the gluon propagator 

D rW = WTk> d, u = g, u -^±^. (60) 

The central question is how to handle this unphysical pole when k + = 0. There are two 
different ways to treat it which employ the Cauchy principal value (PV) and Mandelstam- 
Leibbrandt prescriptions (ML) fl36| : 

PV fc7 ~ 2 {(kn) +i0 + (kn) - io] ' ^ 

ML-L = , (62) 

k + (kn)(kn*)+%0 v ; 

with an arbitrary four- vector n* satisfying n* 2 = 0, nn* = 1 (without loss of generality, we 
can put it equal to p). 

Here, we outline the one-loop renormalization program for the abelian gauge theory with 
PV prescription. In the next section, we show, using a simple example, the difference one 
encounters in dealing with the ML prescription. 

Due to an additional power of the transverse momentum k± in the numerator of the 
density matrix of the gluon propagator, there exist extra UV divergences of the Feynman 
graphs which are absent in the usual isotropic gauges. For our practical aims, we limit our- 
selves to the calculation of the one-loop expressions for the propagators and vertex functions. 
This is sufficient for reconstruction of the equations in LLA using the renormalization group 
invariance. 
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The unrenormalized fermion Green function is given by the expression 

G~ 1 (k) =J£ — m — (63) 

where is a self-energy operator. Calculating the latter to the one- loop accuracy we get 
the following result 

G(k) = (1 - ^U^W—^—U^k), (64) 

Til 



where 



and 



/77 ! 

U ± (k) = 1 - — £ 2 (A;) 7+ - T-(S 2 (A:) - £i) 7+ #, (65) 
fC_|_ /c_^ 

777 I 

tf 2 (fc) = 1 + — £ 2 (A;) 7+ + t-(S 2 (A;) - E x ) ^ 7+ , (66) 
S^^lnA 2 , S 2 (A;) = £lnA 2 /^ ? ^-e? 1 (/,,'-,). (67) 



Here m is a renormalized fermion mass related to the bare quantity by the well-known 
relation 

m = m (l-3Si). (68) 
The functions 6™j 2 ...i used throughout the paper are given by the formula 

_ oo ^-« m n(«^- i +^)- jfc . (69) 

For our purposes it is sufficient to have an explicit form of the simplest functions 

eS(x) = o, (70) 

e° 2 (x)=6(x), (71) 

9^,,,^^^-^-^ (72) 

x 1 - x 2 

since the others can be expressed in their terms via the following relations: 

e° 21 ( Xu x 2 ) = _ ^—e?^!,^), (73) 

Xi - x 2 

e 1 21 (xi,x 2 ) = -J—e?^,^) — ejftxi), (74) 

Xi - X 2 Xi- x 2 

Q° 22 (x 1 ,x 2 ) = - j— ^^ Q\ 1 {x 1 ,x 2 ), (75) 

®lu{xi,x 2 ,x 3 ) = — — — Qu(x 2 ,x 3 ) — — e^^i,^), (76) 

x 1 -x 2 Xi - x 2 

Q\ u (x 1 ,x 2 ,x 3 ) = @u(x 2 ,x 3 ) <~>° u (x 1 ,x 3 ). (77) 

Xi - x 2 x 1 - x 2 
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As we can easily observe, the renormalization constants are not numbers any more but 
matrices acting on the spinor indices of fermion field operators. Moreover, the renormaliza- 
tion constants depend on the fractions ki + . The origin of this can be traced back to the lack 
of the rescaling invariance 

d^v(pk) ^ d^k), (78) 

obeyed by the unregularized propagator. 

An abelian Ward identity leads to the equality of the renormalization constants of the 
gauge boson wave-function and a charge Z 3 = Z g so that the corresponding logarithmic 
dependence on the UV cut-off cancels in the sum of these two contributions in the evolution 
equation and, therefore, we can neglect the fermion loop insertions into the boson line (in 
the QCD case this is no longer true). 

One can easily calculate the vertex function to the same accuracy. The result is 

T p (k u k 2 ) = (1 + S 1 )C/f 1 (A; 1 )^C/ 2 (A; 2 ), (79) 

where 

Q P = l P -{h- m)Q p {k 1 , k 2 ) 1+ - J+Qpih, k 2 )(j£ 2 - m) (80) 

and 

Q p (h, k 2 ) = Es^xh-T+Tp + £ 3 (fc 2 ) 7 p 7+7-, ( 81 ) 

here 

£ 8 (*0 = ^ lnA V dz , (^±e 1 111 (z',z / - z u z' - z 2 ). (82) 

Apart from the graphs we accounted for, there exists an additional UV divergence of the 
virtual Compton scattering amplitude; however, we do not need its explicit expression for 
our practical purposes. This completes the consideration of virtual corrections which cause 
the logarithmic dependence on the UV momentum cut-off of the quantities in question. 



2.7 Difference between the PV and ML prescriptions. 



In the subsequent discussion we will use both the prescriptions for the gluon propagator 
in our practical calculations. The first one (PV) will be used in the momentum space 
PO , EBL |37|, P5| , while the second one in the coordinate space formulation [^, As a by- 
product we verify that both of them do lead to the same result. However, it is worthwhile 
to realize the distinctive features one faces in the computation of the same quantities. 

The most important difference of the second prescription is the presence of the additional 



absorptive part |39| of the vector boson Green function, namely 



Disc 



(k) 



k 2 + i0 



-2m9(k A 



i pu {k)5{k ) - — {k p n v + k u n p )5{k + k_ 



(83) 



The second contribution is of the "ghost" type since it has the wrong sign as compared 
to the conventional one. However, it is not an optional choice but it is an unavoidable 
consequence of equal time canonical quantization [Q. The consequence of this addendum 
can be easily recovered in the calculation of the one-loop evolution kernels. Let us consider, 
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for instance, the coordinate space formalism |33], |38|]. Then, the operator vertices V(Xj) have 



only exponential dependence on the position on the light-cone 



V(Xj) = T v e 



54) 



A simple calculation of the one-loop diagram for the leading twist density with Ty = 7+ 
gives 



7+e 



\Xk-i 



= ^-lnA 2 [ dy 

A 2 2vr Jo 



7+e 



1-2/ 



(85) 



so that the well-defined 1 /^-distribution appears already in the contribution with "real- 
emission". The self-energy insertions into external legs produce a familiar ^-^S(y) term. 
Contrary to this, as we have seen in section [2.6| , the renormalization constants of the field 
operators in the light-cone gauge with the PV prescription turn out to be momentum de- 
pendent, and as a by-product the plus-prescription fulfilling occurs only in the sum of the 
real and virtual corrections. 



2.8 Evolution equations in the abelian gauge theory. 



In this section, we present a pedagogical illustration of the renormalization group mixing 
problem for the redundant basis of unpolarized correlation functions defined by Eqs. (|35[) - 
(|38|), in the framework of the abelian gauge theory. Our aim here is to show the self- 
consistency of the whole approach we have used since the equations derived below satisfy the 
constraint equality given by Eq. (|39|), which are further employed to reduce the overcomplete 
set of correlators to the basis of independent functions. 






(a) 



(b) 



(c) 



Figure 2: One-loop radiative corrections to the two-particle correlators in the abelian gauge 
theory. The fermion propagator crossed with a bar on diagram (c) means the contraction of 
the corresponding line into the point. 

The one-loop Feynman diagrams giving rise to the transition amplitudes of two-particle 
correlation functions into the two- and three-parton ones are shown in Fig. |2| (a,b). The 
last figure (c) on this picture is specific of the vertices having non-quasi-partonic form |5£j] . 
that is for e(x); it displays the addendum due to the contact term that results from the 
cancellation of the propagator adjacent to the quark-gluon and bare vertices. As an output 
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the vertex acquires the three-particle piece. The radiative corrections to the three-parton 
correlators are presented in Fig. |3| (a,b,c). 



(////) 
I it 

(a) 






(d) 



Figure 3: The one-loop renormalization of the three-parton correlation functions. Self- 
energy insertions into external legs are implied. 

A straightforward calculation yields the evolution equations for the spin-independent case 

P 



in the form [38 



M(x) 
e(x) 



i)(x,x') 



+ 
+ 
+ 



(x-P) 



e u n (x,x-(3) 



P 



(86) 



^ / dp[e(P) S^e°n(x,x- (3) + l -5(P - *)} 



- M{/3){2 







(x-P) 



e° u (x,x-p) 



+ xQ\ x (x, x-P) + 26?! (x, x - P) | 



dp'D(P,p'){2 



P 



(x-P) 



e° n (x,x-p) 



+ 



X 



x — P 



e° ni (x,x-p,x-p + p r ) 



+ S(P-x) Jdp''^Ql 11 (P'',P''-P,P''-P') + 2Q° n (x,x-P)^ ] j. 



J7) 



'27r' 



X 



—e(x) - M(x) 



x 



Q° n (x',x' -x) 



X 



-e(x') - M(x') 



x 



+ / dp' (d(x, P') X + *' ) e° ni (x r , x'-x,x'-x + p') 



X 



x'-P 



(x — x') 

-\D(x -x' + P', p') - D(x, x')}Q° u (x', x' - p') 



J dp^D(p, 



(P-x' + x) 



. i ) 1 1 1 ( ^ 3 ^ x ^ x x ~\~ j3 1 



[X — X 



-[£>(/?, x'-x + P)- D{x, x')]Q° n (x, x-P) J- ^D(x, x') 



(88) 



where the dot denotes the derivative with respect to the UV cutoff ' = A 2 d/dA 2 and the 
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plus-prescription is defined by the equation 



P 



(x-P) 



P 



(x-P) 



e° u (x,x-p)-5(p-x) I dp" f e° n (p",p"-p). 



We have used also the equation 



PV / dp 



x 



(x-P) 



Q° u (P,P-x) + Q ( i 1 (x,x-P) 



(P" ~ P) 



0. 



(89) 



(90) 



By exploiting the relation provided by the equation of motion, we can easily verify that 
the RG equations thus constructed are indeed correct and the renormalization program can 
be reduced to the study of logarithmic divergences of the three-parton Z(x, x') and quark 
mass M(x) correlators in perturbation theory. 



2.9 QCD evolution of the twist-3 distributions. 

For the non-abelian gauge theory the equality of the renormalization constants Z g = Z% 
implied above no longer holds; so we should account for the renormalization of the gluon 
wave-function as well as for the renormalization of charge explicitly. For these purposes, to 
complete the renormalization program outlined in the preceding section |2l], we evaluate the 
gluon propagator to the same accuracy. The result can be written in the compact form 

D, v (k) = (l + n*(*)) U„(k)^Z& Um ,{k), (91) 



where 



and 



EW*) = fa ~ l -U add {k) k » n » + Kn » (92) 



U tr (k) = 2— In A 2 \.Ca [ dz ^ , * C e° n (z, z - () - ^ 



4vr 1 "J z(z-()e iiV ^ 3 



Tl add (k) = —\nA 2 C A J dz^^ z(z-()e Qu^^-O (93) 

are the transverse and longitudinal pieces of polarization operator. The renormalized charge 
is given by the well-known "asymptotic freedom" formula 



a , a2 [N f 11^ 
1 + — In A -f - —Cj 



4?r V 3 6 



(94) 



Now we are in a position to adduce the RG equations for the real QCD case. Just giving 
the final result (without intermediate steps) for the chiral-even distributions we then address 
in grater detail to the chiral-odd evolution. 
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2.9.1 Evolution of the chiral-even distributions. 



Let us begin with the transversely polarized structure function g 2 . In the first papers [|I 
on the logarithmic Q 2 -variation of the moments of g 2 a simple evolution has been derived: 



dxx n g 2 (z, Q) 



' <x(Q) 

MQo) 



Jn/Po 



dxx n g 2 (z,Q Q ), 



(95) 



and the anomalous dimensions have been found by calculating the radiative corrections to 
the operator 



O 



afli/12-..^r, 



i n A S ipwsD^D 



1-1-2 



D»J> 



(96) 



over the free quark states. As we have seen above, the conclusion of these works is erroneous 
since (j96|) mixes with other operators of the same twist and quantum numbers due to renor- 
malization. The most significant departure from the naive expectation is that the number 



of operators involved in the RG evolution turns out to be increasing with the moment |3C 



and as a result it is impossible to write the equation of the DGLAP type which manages 
corresponding scale dependence for the g 2 . Moreover, the eigenvalues of the anomalous di- 
mension matrix are not known analytically. However, an way out has been found |32| in 
the two important limits: N c — > oo and x —* 1, where evolution reduces to the DGLAP 
equation (1951), though the anomalous dimension turns out to be different from that found in 



40|. The combining use of these asymptotics does provide an excellent approximation [B8 



to the complete result. 

To simplify the discussion, we neglect the quark-mass operator. The nonlocal string 
operators ^Sp introduced previously are related to each other by charge conjugation (so are 
the corresponding kernels which govern their evolution) and do not mix in the course of 
renormalization; therefore, we give only the result^ for ~S a 



S a (ii,X) 



a 



, dy [ dz{ Ca 
2tiJo y Jo 1 



{N(y,z)] + - 7 -5(y)5(z) 



S a {jj,y, X- jjt,z) 



2z+[N(y, z)] 
C A 



+ 

+ (CV- 



S(y)S{z) 



S a (/i - Xz, Xy) 



y5(z) S a (-[My, X-fxy) -2z <S CT (//- Xz,-Xy) 



+ i K (y, z ))+ S a {jJ,z + Xz, Xy + fiy) 



(97) 



To condense the notation we have used the dot as short-hand for the logarithmic derivative 
with respect to the renormalization scale ' = fi^d/d^ and ~S p (fi,X) = ~S p (fi, 0, A). The 
standard plus-prescription fulfilling is 

[N(y, z)} + = N(y, z) - 5(y)6(z) f dy' f dz'N(y', z'), N(y, z) = 5(y - z) V - + S(z)^ 

Jo Jo y y 

-l 



[K(y, z)]+ = K(y, z) - 5{y)5{z) f dy' T dz'K(y', z'), K(y, z) = 1 + 6(y)~ + 8{z) V -. 

Jo Jo z y 



(98) 



This discussion follows Ref. [B3| 
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To proceed further we construct a C-even operator from the Shuryak-Vainshtein ones 



y a (\, f j l ) = + sf > (x,fjL) + -s f, (ji,x) 



(99) 



and define a new distribution function as Fourier transform with respect to the variable A 
only, so that 

y(x,u)S$ = I [ ^e iXx (h\y a (uX, -uX) + (A -> -X)\h) (100) 



2 J 2tt 

depends on the effective momentum fraction x and the gluon position on the light cone u. 

Skipping the details (to which we address in the following discussion of the chiral-odd 
distributions), we just note that the genuine twist-3 part of g^ which can be expressed via 
the integral of the three-parton correlator 



92[X) 



with 



~92{x) + j —gi(x) 
Jx y 



xg2(x) = — — / duuy(x } u 
ax Jo 



(101) 



(102) 



satisfies the DGLAP evolution equation in the large- iV c limit (i.e. neglecting the terms 
0{l/N 2 )) 



[xg 2 {x)\ 



a 



4nJ x y 99 



x 



[y92{y)}, 



(103) 



with the splitting function ||32| , [33 



P 99 (y) = n c 



! 


"2" 




y. 



2 - y + 



(104) 



In the x — > 1 region the evolution equation remains of the DGLAP type even for the jj- 
suppressed contribution. The combining use of the Eq. (|104[) and an additional piece 



AP„(y) = - K 



(105) 



for d?(l/iV c )-terms with x — 1 1 yields a good approximation for the exact evolution. Ob- 
viously, the function gz(x) obeys the same evolution equation as g~i. Actually, the splitting 
function (|104[) has been exploited in Ref. |13|] to rescale the bag model predictions to values 
of Q 2 of the real experiment. It has been shown there that g% (x) and g%{x) enter into 
g2{x) on equal footing at the model scale //bag- Moreover, it has been figured out that this 
situation does not changed at higher Q 2 and g2{x) remains an important ingredient of g2{x). 



2.9.2 Evolution of the chiral-odd distributions^. 

Turning to the case of the chiral-odd distributions, we have to note that in the leading 
logarithmic approximation the evolution equations that govern the Q 2 -dependence of the 

5 The following discussion mimics the corresponding sections from Ref. J3S| 
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three-particle correlation functions are the same, discarding the mixing with the quark mass 
operator. Therefore, we omit the "tilde" sign in what follows. 

In the light-cone fraction representation we get for the correlation function D(x,x') 



D(x,x') 



a 

'2lT 



-a 



(x — x') 



x'M(x)G n (x', x'-x)± xM(x')Q n (x, x - x 1 ) 
C 



+ j dp\C F D{p, x')^&° n (^ x-x') + ^ ([D((3, x') - D(x, x'^-^-^Q^x, x-P) 



+ [D(p + x',x') - D(x,x')\- 



(x — x 



x — x' — P) 



x j x x ^3 ^ 



+ 



(P+x- x') 



X 1 



(d(/3, x') X _ e^i, x - P) + D{fi + x', x')e° n (x -x',x-x'- p)X\ 



c f 



c f 



D(/3,x 



ls ( P + x- x') 
(V — x) 



Q® u (x, x — x', x — x' + (3) 



[D(P, x'-x + P)- D(x, x')]^—-G° u (x, x-P) 



+ 



J dp' (c F D(x, P'^e^x', x'-x) + ^f ([D(x, P') - D(x, x')} -^^ e u (x', x' - P') 



- [D{x, P' + x)- D(x, a /)] 7 J^-j— e° n (x' -x,x'-x- P') 



(x'-x- P') 



+ 



(p' + x'-x) 



X 



D(x, P')^- — rQUx', x 1 - p') + D(x, P' + x)e u (x' -x,x'-x- p') 
[x — x') 



C 



+ [ c f-^t) [D(x,p') 



?/ . [P' + x'-x)^ 



[x — X) 



11 1 i^x ^ x x j x x ~\~ j3 ^ 



+ [D{x-x' + P' ) P')-D{x 1 x')}-^ i Q\ 1 {x' ) x' -P')^\ -^C F D(x,x') 



P' 

and for the mass-dependent correlation function we have 



M(x) 



-Ci 



a 

'27T 



dpM(p) { 2 



P 



(x-P) 



e° n {x,x-p) 



+ ^&°n(x,x-P) 



(106) 



(107) 



where we have used the standard plus-prescription fulfilling / dx [...]+ = (for a definition 
see Eq. (p9|)). Throughout the paper the plus and minus signs in the mass-operator term 
correspond to the functions D (for e) and D (for hi), respectively. 

For the string operators we obtain the following compact RG equation: 



Z(\,y) = m dy /_" dz{C F y 2 5(z) M 1 (A - fiy) ± M 1 (Xy - fi) 



+ 



2z+[N{y,z)) + -^6{y)S(z) 



[C F -^ 



[L(y,z)] + --5(y)5(z 



2z [Z(-Xy, n-\z) + Z(X - fiz, ~ny)] 



[Z(Xy, /I- Xz) + Z{X - fMZ, ny)\ 
Z(Xz + \iz, \iy + Xy) 



(108) 
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with the function N(y, z) given by Eq. (|98|) and 



z y 



The equations written so far should be supplemented by the following 

-2-y-y 2 )M\Xy), 



a 

— ( 

2tt 



a. 





~2 








V. 




"2" 








y. 



2 + -5{y)\h x {Xy). 



(109) 

(110) 
(111) 



The last one, when transformed to the momentum space using the formulae of the next 



section, coincides with the result obtained in Ref. |L6 
2.10 Local anomalous dimensions. 

Now we are able to pass from the evolution equations for correlators to the equations for their 
moments and find, in this way, the anomalous dimension matrix for local twist-3 quark-gluon 
operators. 

We define the moments in following way: 

F n = J dxx n F{x) for any two-particle correlator, 

Z l n = J dxdx , x n - l xf l - 1 Z(x,x'). (112) 

In the language of operator product expansion these equalities specify the expansion of 
nonlocal string operators in towers of the local ones, namely 



Z, 



2—1 Qn~l 

— 7 Z(A,//)| A=M=0 



dn 1 - 1 dX r 



(iD-i 



\n—l 



M n 



M(X)\ 



A=0 



(113) 



The inverse transformations to the nonlocal representation are given by 

co m \n oo 

Now it is a simple task to derive the algebraic equations for the mixing of local operators 
under the change of the renormalization scale from the evolution equations ( |106|) - (|110|) . They 
are 

a 

~Z MM 

In 



{ J n M n 



Zl = — 
2tt 



ZMl n n _ l+1 ± ZMl ^ M n + J2 

zzl^kZ, 

k=l 



n ryk I 
n ( ' 



(115) 
(116) 
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where the anomalous dimensions are given by the compact expressions 



mm7 
zmIi 



— Cf (S n + S n +2) > 

2C F 
Z(/ + l)(/ + 2)' 

\c F 5(l - k) + 9a 



;ii7) 
;ii8) 



9(1 



(k + l)(fc + 2) 
(Z-fc)(Z + !)(/ + 2) 



5(1 -k) [Sk-i + S k+2 ] 



+ 



+ 6(1 -k) 



9(1 - k - 1) 
2(-l) fc 



2(-l) fc Cf 



1 



ii— k 



k(k + l)(k + 2 
Here, we have used the following step functions: 



l(l + l)(l + 2) (l-k) C 1 - 1 
k — ► n — + l\ 
/ -> n-Z + 1 J' 



5, 



+ 



;ii9) 



1, i > j 
0, i < j 



6{i-j) 



1, i=j 
0, 



(120) 



The 



as well as the convention S n = Y2-i f and the binomial coefficients C" 1 = ,, nl — rr. 

" *—iK— i n m!(n— m)! 

results of this section have been independently derived in Ref. [[Il]] using the standard 
approach based on the local operator product expansion. 



2.11 Relating the evolution kernels in the light-cone position and 
light-cone fraction representations. 

Until recently the relation between different formulation of the evolution equations in the 
light-cone fraction [|30], and light-cone position [27, 0] representations has been obscure 



and has been thought to be difficult to realize p3| . However, having at hand the evolution 
equations in different representations for the same quantities, we are able to fill this gap and 
relate the kernels in both the cases [j38| . Such a bridge can be easily established using the 
Fourier transformation for the parton distribution functions given by Eqs. (|3[) and (|TT|). 

To start with, we address ourselves to a simpler case of the two-particle correlation 
functions T. The evolution equation in the light-cone position space is of the following 
generic form: 

^(A)= I' dytC(y)F(\y), (121) 
Jo 

where JC(y) is the corresponding evolution kernel. By exploiting the definitions @ we can 
recast the Fourier transform on the language of two-particle evolution kernels. In this way, 
we find the direct transformation 

K(x,P)= [ 1 dyJC(y)5(x-y^. (122) 

With the help of the general formula 



o 



Q ' dyf(y)5(x -y0) = f B° n (x, x-(3). (123) 
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we observe that the RG equations for two-parton correlators derived in the previous section 
are the same indeed. The inverse transformation can be done 



*^ e -o*W K ( $ = ( X dyK{y)5{n-y\) 



with the following result 



2n 



dyf{y)S(fji-yX). 



(124) 



(125) 



The transformation for three-particle correlators is a little bit more involved. The general 
form of the evolution equation for the light-cone string operator Z(\,fi) reads 



dy / dzK(y, z)Z(r) U \ + i] 12 fi, r] 21 X + r} 22 \i), 



(126) 



where 77^ are linear functions of the variables y, z. In the momentum fraction representation 
the evolution equation looks like 



Z(x, x') 



d/3d/3'K(x,x',/3,/3')Z(/3,/3'). 



(127) 



Specifying the particular form of the functions rjij, we display below an example for the 
conversion supplied with a general formula suitable for all practical cases of interest. 
For the C^/2 part of the evolution equation, the Fourier transformation gives 



:i28i 



K(x, x', (3, (3') = 8{(3' — x') I dy I dzJC(y, z)5(x — x'z — (3y). 
The particular contribution is 



IC(y,z) = z ^ K(x, x', (3,(3')= 5((3' — x') | — j^-Ex(x,x — x',x — (3) 



x 



+ ^-E 2 (x,x - x',x- 

X J 



Here, we have used 



and the following general result: 



E n (x, x-x',x-(3)= dyy n e° n ((x - (3) + y(3, (x - (3) - y(x' - (3)) 



(129) 



(130) 



1 

n 



1 - 



' (3-x 
B-x> 



e° u (x,x-x') + -^ 
nx' 



13 — x 
0-x> 



'(3-x 



e° u (x,x-p). 



The complete list of transformations can be found in Ref. |38 |. Using these results, we 



can easily verify that the evolution equations given by Eqs. 



and 



agree with each 



other. It should be noted that it is sufficient to have at hand Eqs. ( |123|) and ( |130| ) to perform 
the conversion from one representation to another. 
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2.12 Generalized DGLAP equations for the three-parton correla- 
tors. 

As we have seen above, the evolution equations in the momentum space ( |106| ) turn out to be 
very complicated, while Eq. ( |108|) being compact is not suitable for an analysis since only its 
Fourier transform is related to the physical observables. Therefore, to obtain a simple and 
manageable equation which can be attempted to be diagonalized, we are forced to proceed 
further in the same line as suggested in section |2.9.1 for the chiral-even structure function. 
For this purpose we define a new function, Fourier-transformed with respect to the A variable 
only: 



z(x, u) 



~J^(h\Z(u\,-u\) 



± (u 



u)\h), 



which is even under charge conjugation and depends on the variables x and u. The latter 
has the meaning of the relative position of the gluon field on the light cone. For < u < 1 
the gluon field lies between the two quark fields. Because of the support property \x\ < 
max(l, \2u — 1|), the variable x is then restricted to \x\ < 1 and can be interpreted as an 
effective momentum fraction. 

The evolution equation for z(x,u) can be derived in a straightforward way from the RG 
equation ( p.08| ) for the nonlocal string operator Z. It can be presented in the form of a 
generalized DGLAP-type equation in the mixed representation 6 : 



z(x,u) = ^- J y J dv^P zz (y,u,v)z(^-,vj + P zm (y,u,v)m (^j 
dy 



m(x) 



a: 
2^ 



X 



(132) 
(133) 



Here, m(x) = xM(x) and the integration region is determined by both the support of z(x, u) 



and the kernels 
' C F 



x, u, v) 

Ca 
2 



6>i(x, u, v) [L(x, u, v)] + — 2 (x, u, v)M(x, u, v) — -5(u — v)5(x) 



M(x, u, v) + [N(x, u, v)] + — -8(u — v)5(x) 











\v - 





p. 



X 



Cpx 6(x)6(x)— [8(v — it — xu) ± 5(v 

V 



u 



xu)} , 



p„ 



x 



C f 



(134) 
(135) 

(136) 



where the auxiliary functions are defined by: 

0i(x,u,v) = 9(x)9(u — xv)0{u — xv), 

6 This particular representation of the evolution equation for three-particles distributions has first been 
given in the second paper of Ref. |33| . 
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2 (x,u,v) 
3 (x,u,v) 
L(x, u, v) 
M(x, u, v) 
N(x, u, v) 



u J \ u J \ u 



UJ \ u 
u 2 



v(v — u) 
2x(l — xv 



5(u — xv ), 



M 3 



ve(u) 



— o(x) H o[u — xv) 

U V 



(137) 



u[v — u) 

The plus-prescription for the arbitrary function A is defined by the equation 

Oi(x,u,v)[A(x,u,v)] + 

= 0{(x, m, v)A(x, u, v) — S(x)5(u — v) J dx' J dv' 0i(x',u,v')A(x',u,v'). (138) 
Note that, due to the evolution, the variable u is no longer restricted to the region < u < 1. 



2.13 Eigenvalues and eigenf unctions. 

Going further, for simplicity, we restrict ourselves to the homogeneous case, i.e. we discard 
the quark-mass operator, which is certainly a justified assumption for the light u- and d-quark 



species. The diagonalization of the evolution equation ( |132|) can be achieved by introducing 
the Mellin transforms 

„n-l. 



Z ( It 



dxx n z(x,u), 



(139) 



where n is the complex angular momentum. Operators with different n do not mix with 
each other and satisfy the following equation: 



z n (u) 



a 



With the kernel given by 



P n zz M 



+^0 3 (u,v) 



C, 



2 J 



_ } , drPUa.r) 



e x (u, v) [L n (u, v)]+ - 2 (u, v)M?(u, v) - -5(u - v) 



(140) 



7 



M 2 n (w, v) + [N n (u } v)} + - -5{u - v 



+ 



u — > u 



(141) 



where the auxiliary functions read 

G 1 (u,v) = 9(v-u), 2 (u,v) =G(-v)0(l-vu), G 3 (u,v) = 9(v)9(v - u). (142) 



L n (u,v) 
M?(u,v) 



e(v) 



v — u \v 
2 f 1 



n+l 



u 3 L n + 1 



,71+1 



u 



n+l 



n + 2 



,n+2 



U 



n+2 
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The plus-prescription is defined as 



Oi{u,v)[A n {u,v)} + = Oi{u,v)A n {u,v) - 5{u - v) J dv'Oi{u,v')A n {u,v'). (144) 
Here, we can check our calculation once more since in multicolour limit, exact Eq. ( |140| ) is 



reduced to the approximated equation of Ref. [44 



To obtain the solution of the evolution equation ( |140| ) , we choose n as a positive integer. 
In this follows from the definition ( |131|) of z n [u) the n-th moment is actually given 

by the following linear combination of local operators Z l n (see Eq. (|112|) ): 



n 

-E c ti« B " , « / " 1 ^i. ( 145 ) 



Z (U 



1=1 



so that z n (u) is a polynomial of degree n — 1 in u. Thus, the kernel P zz (u,v) possesses n 
polynomial eigenfunctions ef(v) 

dvP2 z (u,v)e?(v) = -\fe?(u), Z = l,...,n, (146) 

where — Xf denotes the eigenvalues. The solution we are interested in is given in terms 
of the eigenvalues and eigenfunctions of the anomalous-dimension matrix Z zl l n of the local 
operators Z l n . The eigenvalue problem we have attacked has no analytical solution; however, 
the diagonalization can be done numerically for a moderately large orbital momentum n, 
e.g. n < 100, which is quite sufficient for practical purposes. These eigenfunctions can be 
constructed by diagonalization 

/n 
dvP^v^v"- 1 = Y.Cn-\z Z llu n - l u l -\ (147) 
1=1 



where the anomalous- dimension matrix zz ^ k of the local operators is given by Eq. ( |119| ). 
Actually, this is a purely algebraic task and we find 

n 

e n k (u) =Edrt i_ X> with {{E n Y l zzl n E n } ki = -X n k 5(k - I), (148) 
i=i 

where 5{k — I) is a Kronecker symbol defined by Eq. ( |120| ). The spectrum of the eigenvalues 
A" up to n = 50 is shown in Fig. |4j. The solution for the moments z n (u) (in the massless 
case) is then expressed in terms of the eigenfunctions and eigenvalues we have found p8| 



z n (u,Q 2 ) =±c?(Ql)e?(u)e W l [ -Q j^A?|. (149) 

The coefficients cf {Qq) at the reference momentum squared Q$ have to be determined from 
the non-perturbative input. 
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Figure 4: The spectrum of the eigenvalues Xf for the evolution kernel P£ z defined in (|141 ) 



Obviously, in order to find the evolution of the higher n-moments of the structure func- 
tions the whole information about the relative size of the reduced matrix elements of the 
(I — 1, ... ,n) local quark-gluon operators is needed. At the moments, in the lack of complete 
understanding of the yet unclear confinement mechanism this problem is not accessible by 
nonperturbative methods presently available. However, we should note that in the nearest 
future it will not be possible to distinguish experimentally between the terms with different 
anomalous dimensions even for the transverse structure function g2{x) to say nothing of 
h L (x) (not to mention e(x)). 

Eq. ( |149|) can be rewritten directly for the structure functions entering into the physical 
cross sections. To this end, we have to dispose the relations similar to the ones given by Eqs. 
( 39 ) and fl5"6|) transformed to the mixed representation ( [L31] ). Namely, we have 



e(x) = — / duz(x,u), 

2 dx Jo 

1 d f 1 

h L (x) = — / du (1 — 2u)z(x, u), 

2 dx Jo 

where we introduce for convenience a new function hi(x), so that hi(x) reads: 

-'- d r 1 dyx 2 z 



h L (x) 



x 



dy-zh^y) -x 
y A dx 



y y 



My). 



(150) 
(151) 

(152) 



The last term on the RHS coincides with the twist-3 part hi,. 

Thus we arrive to an explicitly diagonalized form for the structure functions f(x) (sche- 
matically): 



[f{Q 2 )]n 



l duW{u)z n (u,Q*) =£rtit(E? k r 1 (h\Z k n (Q )\h) 



(153) 



1=1 



k=l 



where W^(u) is a weight function in Eq. ( |150| ) and ( 151 ) and J l n is the overlap integral 

J l n = [ duW f (u)e?(u). (154) 
Jo 
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2.14 Solution of the evolution equations in multicolour QCD. 



The complicated form of the evolution ( |153j ) of the twist-3 distributions compels one to look 
for the simple approximated solution which could work with reasonable accuracy. The reso- 
lution of this problem is based on the observation that in the large- N c limit only the planar 
diagrams (Fig. [| (a,d)) survive and the kernel P zz (u, v) has two known dual eigenfunctions: 1 
and l-2u g fS|, so that J* du ef (u) = 5 n +0(l/N c ) and /„* du (l-2u)e]*(it) = 5 l2 +0(l/N c ), 
where I — 1,2 correspond to the lowest two eigenvalues of the spectrum shown in Fig. |], 
thus in the sum in Eq. (|153|) only the lowest two terms survive. Then, a straightforward 
calculation gives the following DGLAP evolution kernels: 



du 



1 

l-2u 
l-2v 



P 2Z {x,u,v) = N c 9{x)d(x) 



_L 1-2 _ 5 
~ 2 4 



3 2 
2 



5{x) 



-Ax) 



(155) 



As we have observed previously in section |2.9.1| in the context of the chiral-even distribution 
g2{x) f3"2]| , similar equations hold true also for the -^--suppressed terms in the x — > 1 limit 
for flavour non-singlet twist-3 evolution kernels. In the present chiral-odd case, we find 



o 



{ l-2v ) 



^9(x)9(x) 



+ I5(x) + <D(x< 

+ m(x) + o(x 



+ N C 



(156) 



where the N c ■ ■ ■ symbolize the x — > 1 limit of Eq. ( |155|) . 

The eigenfunctions we have obtained coincide precisely with the coefficients W" (u) that 
appear in the decomposition of e(x, Q 2 ) and hi,{x, Q 2 ) in terms of three-particle correlation 
functions. 

From the observations we have made above, it follows that in the large-iV c as well as 
in the large-x limit the twist-3 distributions satisfy the DGLAP (ladder-type) evolution 
equations that hold for the twist- 2 operators. By combining the large-iV c evolution with 
the large-x result for the ^--suppressed terms, we can improve the accuracy of multicolour 
approximation within a factor 5-10 and reach the precision of a few per cent as compared 
with the evolution predicted using an exact equation ( |140|) but supplied with a model for the 
light-cone position distribution of gluons between the quark fields [38|. Thus, the functions 
e(x, Q 2 ) and Jil(x,Q 2 ) obey the following improved evolution equations: 



(157) 



with 



P ee {y) = 2C F 
P- hh (y) = 2C F 



C A 



3C 



A 



c 



o C A j 5(y) + O (f/N c ) 

7C F _ AC A ' 
~6 3 



5(y) + O (y°/N c 



158) 



Note that h F {x) fulfills the same evolution equation as hi(x). The simplest way to verify 
this is to make the Mellin transform of the corresponding evolution equations. 
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2.15 Remarks on the momentum space formalism. 



Let us add a few remarks on the momentum space formulation. As we have seen above the 
solution of the evolution equations in the asymptotic regimes is the most straightforward in 
the light-cone position representation. It is by no means trivial to observe the appearance 
of the DGLAP equations in momentum fraction representation. However, we know that the 
asymptotic solution, in coordinate space, is given by the convolution of the three-particle 
correlation function with the same weight function that enters in the decomposition of the 
two-parton correlators at tree level. With this in mind, we are able to check that the integrals 



ex 



h L (x) 



x 



d(3'D(x,(3'), 
1 df3 f dp' 



2 + (P~P') 



d_ _d_ 



D(P,P% 



(159) 
(160) 



taken from Eqs. ( |3"§D and 
DGLAP equations, namely 



P'~P 

neglecting quark-mass as well as twist-2 effects, satisfy the 



e(x) 



(x-P) 



e° u {x,x-p) 



~^N C J d/5e(/3){2 



(161) 



h L (x) 



a 



—N c I dph L (P)\2 



P 



(x-P) 



+ 



x 



2-3-1 e° n (x,x-P) 



e n {x,x-p) 



-S(P-x) 



(162) 



The corresponding anomalous dimensions are 



[hi 



a 

Air' 



-2^(n + 2) - 2 lE + 



n + 2 



[e]n, 
[h L }n- 



(163) 
(164) 



Which are exactly the anomalous dimensions 7^ found in Ref. [44] for e and hL, respectively, 
with the replacement n — > j — 1 Q. 

Concluding this section, we have to note that the simple DGLAP equation ( |162| ) has 
been used, recently, for prediction of the size of the structure function h L (x) |L4| at high Q 2 
starting from its value in the low energy point found in the framework of the MIT bag model. 
The main conclusion of this work is that the twist-3 contribution to hi(x) is significantly 
reduced in the course of the evolution in contrast to the corresponding situation in the case 
of g2{x) structure function discussed above. This observation is the direct consequence of the 
larger anomalous dimensions for hi{x) at low values of n-spins as compared with g2{x). This 
means that it will be extremely difficult to extract hi(x) at large momentum transferred. 

7 The difference in the anomalous dimensions is due to an extra power of the momentum fraction x 
included in the definition of the twist-3 correlation functions. 
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However, if the latter will be sizable at high Q 2 in future experiments it will indicate that 
the naive bag model predictions could not be trusted for the calculations of the quark-gluon 
correlations presented in hi(x) (as well as in the other twist-3 distributions). 



3 Time-like cut vertices. 

3.1 Time-like processes and the Gribov-Lipatov reciprocity. 

The deep inelastic scattering of lepton beam on the hadron target has proved to be the most 
effective experimental tool for studying the dynamics of hadron reactions on the parton level 
which has a firm basis in the quantum field theory provided by the light-cone OPE. As we 
have seen in the preceding sections it makes possible the investigation of the logarithmic 
violation of the Bjorken scaling as well as the power suppressed contributions responsible for 
polarized phenomena. However, we mainly use an equivalent approach for the analysis of 
the corresponding quantities which is based on the factorization theorems and the evolution 
equations since the latter can be applicable in the situations when the OPE is no longer valid. 
These are the inclusive production of the hadron in the e + e~-annihilation, semi-inclusive deep 
inelastic scattering, Drell-Yan lepton pair production et ctr. 

There is continuous interest in the inclusive production of hadrons in hard reactions. 
These processes involve a quark fragmentation function to describe the hadron creation from 
the underlying hard parton scattering. However, they differ considerably from the DIS as 
the short distance expansion could not be employed although the given processes go near the 
light cone. The theoretical basis for strict analyses of the above phenomena is realized by the 
generalization of the OPE to the time-like region in terms of Mueller's £-space cut vertices 
As we have mentioned in the introduction an essential departure from the DIS is that 



the moments of the fragmentation functions are essentially nonlocal in the coordinate space. 
However, this approach has all attractive features of the OPE as it provides a consistent 



framework to account for the higher twist effects [28] as well as it allows to sum up the UV 



logarithms j|o] by using the powerful methods of the renormalization group. 

The semi-inclusive hadron production from a quark fragmentation is described in QCD 
by the specific nonperturbative correlation functions of quark and gluon field operators over 
the hadron states which can be identified with £-space cut vertices. While the behaviour of 
the latters with respect to the fraction of the parton momentum carried by the hadron is 
determined by the nonperturbative strong interaction dynamics, the large Q 2 -scale depen- 
dence is governed by perturbation theory only. Since the cross section we are interested in 
cannot be related to the imaginary part of some T-product of currents, therefore, we must 
deal with particular discontinuities of the Feynman diagrams from the very beginning. 

Inasmuch as the twist-3 fragmentation functions enter into several cross sections on the 
same footing as the distributions, their scale dependence is of great interest. Apart from 
significance for phenomenology, it is important for theoretical reasons: while it is know 
that in the leading order of the coupling constant the splitting functions for the twist-2 
fragmentation functions can be found from the corresponding space-like quantities via the 
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Gribov-Lipatov reciprocity relation : 

PSL {x) = pTL Q\ ^ (165) 

(note that both quantities are defined in the physical regions of the corresponding channels) 
no such equality is known for higher twists. 

We begin our discussion with a study of twist-3 nonpolarized chiral-odd (NCO) fracture 
functions. This is the simplest case with respect to the number of correlation functions in- 
volved in mixing under renormalization group evolution. From the phenomenological point 
of view they appear, for example, in the cross section for semi-inclusive hadron (H) produc- 
tion in the process of measuring the nucleon's (h) transversity distribution hi(x) from deep 
inelastic scattering [ |H|j : 

d*Aa 4oL r A V\ n ( r\ 



em 



dxdyd(l/()d(f) Q 2 

+ cos sin ^(k- 1)(1 ~y) (g t (x,0-G 1 (x,0 (* ~ §))]• ( 166 ) 

Here k — 1 + Ax 2 M 2 /Q 2 , y = 1 — E'/E and the cross section is expressed in a frame where 
the lepton beam with energy E defines the z-axis and the x — z-plane contains the nucleon 
polarization vector, which has the polar angle \ and the scattered electron E' has the polar 
angles 9, 0. The functions G\ and Gt are expressed in terms of the product of the familiar 
distribution and fragmentation functions in the following way 

G i(^C) = l^Q-gl^iO, (167) 



G T (x,0 = \Y,Q 2 i 



g l T (x)r(C) + -h\(x)T(C) 

X 



(168) 



All of them have the expressions in QCD in terms of the light-cone Fourier transformation 
of correlation functions of fundamental quark and gluon fields over specific hadron states 



10| , [17[ . Some of the definitions were given already by Eqs. (0), while the others are 
Sihi{x) = y^e^(h\mK+^(^)\h), 

HO = ^J^(OM(Xn)\H,X)(H,X\m\0)- (169) 

Note that we have used slightly different definition of the function h\(x) which has another 
dimension in mass units as compared to previous one. This is done in order to make the 
product hi{x)X{C) dimensionless. Of course, they can be made dimensionless separately 
by introducing certain characteristic scale rn 2 ^ (which can be set equal to the mass M of 
the appropriate hadron) into the definition of the corresponding correlators. 1 is given by 
equation (|170|) . Note, that the physical regions are different for distribution and fragmenta- 
tion functions: < x < 1 and 1 < Q < oo, respectively. 

In the following sections we address ourselves to the problem of building of the master 
equation for the function X. This problem is of the same complexity as for corresponding 
quantities in the space-like domain since we face the mixing with other cut vertices of the 
same twist and quantum numbers in the course of the renormalization group evolution. 
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3.2 Recombination functions and their support properties. 

For our purposes it is much more suitable to deal with correlation functions listed below 
which are the generalization of the formulae given in section pT3| to the fragmentation region. 




Figure 5: Graphical representation of the three-parton recombination functions. 
Diagrammatically, the three-parton correlators are presented in Fig. |5|. 

J(C) = \j^ X( (OmXn)\H,X)(H,X\m\0), (170) 
M(() = ±f ^e^(Q\my + i>(\n)\H,X){H,X\m\0), (^1) 
ZfHCO = ^J ^^e^-^\0|g 7+7 ^(An)| J ff,X)(F ) X|V;(0)^(H|0), (172) 



4 2 \CC) = }-J ^ge^'^(0|g 7p V^(H^(0)|^X)(i7,X|^(An)|0), (173) 



4 1] (C,C) = ^ / ^^e^'-^(0|^(0)g 7+7 X/in)|^X)(i/,X|^(An)|0), (174) 

4 2 \CX) = ^J ^ge^-^\0|^(An)|F,X)(if,X|^Mg 7 ^(0)|0). (175) 

The summation over X is implicit and covers all possible hadronic final states populated 
by the quark fragmentation. Again we use the light-cone gauge B + = 0, otherwise a link 
factor should be inserted in between the quark fields to maintain the gauge invariance. The 
quantities determined by these equations form the closed set under renormalization, however, 
they are not independent since there is relation between them due to equation of motion for 
the Heisenberg fermion field operator 

1(C) - M(0 - J dC%(C,0 = 0. (176) 
Here and in the following discussion we introduce the convention 

ZjiCO = \ \zf\dX) +-zf(C,C0] • (17T) 

While the former two functions X and Ai can be made explicitly gauge invariant by inserting 
the P-ordered exponential (which is unity in the gauge we have chosen) between the quark 
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fields, the latter can be written in the gauge invariant manner by introducing the following 
objects: 



ft! (C',C) = C%(C',C), ^2 (CO = C^(C,0- 



'1781 



Taking into account Eq. (f|) it is easy to verify, in the same way as before, that they 
are indeed expressed in terms of correlators involving the gluon field strength tensor. The 
functions Z^ and Z^ are related by complex conjugation 



^(C, ' = z?\c, CO, CO ' = z^iC, 



r(2) 



'(2)/ 



(179) 



Their support properties can be found by applying Jaffe's recipe |24 ]. It has been shown that 
the field operators entering the definition of the correlation functions can be placed in the 
arbitrary order on the light cone with appropriate sign change according to their statistics. 
Then taking the particular ordering and saturating the correlation function by the complete 
set of the physical states we immediately obtain (for definiteness, we consider the function 

4 l) ) 



4 1] (C,o 



i 



l-Cx)5(C-Cy)(0\m,X)(H,X\^\Y)(Y\B ± \0) 



- £ 5(C - 1 - Cx)S(C -('- C Y )(0\m,X)(H,X\B ± \Y)(Y^\0), (180) 

4< 5 X,Y 



with (x,(y > and we omit the unessential Dirac matrix structure of the vertex. From 
these equations the restrictions emerge on the allowed values of the momentum fractions: 
1 < C < oo, < < (. By analogy one can easily derive similar support properties for 
other functions. 



3.3 Feynman rules for the discontinuities of the diagrams. Kel- 
dysh diagram technique. 

It is well known that the time-like cross section could not be related to the imaginary part of 
any T-product of the currents. Contrary, it is given by the particular absorptive parts of the 
ladder-type diagrams. In the physical region of the annihilation channel these graphs possess 
additional discontinuities, which are not relevant for our purposes, since we are restricted 
to the cuts that separate the possible hadrons in the final state. To be able to extract the 
imaginary part we are interested in we have to label in some way the field operators in the 
amplitudes to the right and to the left of the cut. This can be suitably done with the help 
of the Keldysh diagram technique [|46[| . It allows to recast the program for the calculation 
of the particular discontinuities of the Feynman diagrams to the operator-like language^]. 
Consider, for instance, certain S'-matrix element which is given by the following functional 
integral 

M = J D<P (0!0 2 ...0„) exp (i J dzL{<f>)\ , (181) 
8 In this discussion we will follow Ref. |2§| 
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where = ip,ip, B. The cross section a = MM* of the process is 



MM* 



J D + (f)D 



exp ji J dz + L((p) - i J dz ~X(» j • (182) 



Here the "plus" and "minus" superscripts label the fields from the direct and conjugated 
amplitudes. Following the original works ||6| one can accept that they are the components of 



a unique operator <P(z,t) composed from the time- and anti-time-ordered fields, i.e. <p + and 
respectively. Now, the Green function for the "big" field is a 2 x 2-matrix constructed 
from the usual Feynman propagator, its conjugated analogue and its discontinuity via the 
Cutkosky rules for the lines connecting the direct and final amplitudes. Thus, the radiative 
corrections to the bare cut vertex can be calculated then by using the conventional Feynman 
rules with the following modifications: 

• All propagators and vertices on the RHS of the cut are Hermitian conjugated to that 
on the LHS. 

• Every time crossing the cut the propagator l/(k 2 — m 2 + iO) has to be replaced by 
-2m5(k 2 -m 2 ). 

• For each propagator crossing the cut there is a ^-function specifying that the energy 
flow from the LHS to the RHS is positive. (In the infinite momentum frame this is the 
plus component of the four- momentum.) 

These statements complete the rules to handle the cut vertices. They can be summarized in 
the Fig. g, where V^ up = (k x - k 2 )^g U p + (k 2 - h) u g w + (k 3 - ki) p g^ u . 



3.4 Abelian evolution. 

In this section accepting the diagram technique derived in the preceding section we show 
our machinery on a simple example of abelian evolution and generalize it afterwards to the 
Yang-Mills theory. As in |2l| we start with overcomplete set of the cut vertices defined by 
Eqs. (|170|) - (|175|) and disregard for a moment the relation between them. Then Eq. (|176|) 
verifies that the evolution equations thus obtained are indeed correct. We have to note that 
since the observed particle H is always in the final state some cuts of Feynman diagrams 
are not allowed and, therefore, we could not obtain the evolution kernel for the cut vertex 
taking naively the discontinuity of uncut graph as we are restricted over the limited set of 
the cuts. 



3.4.1 Sample calculation of the evolution kernels. 

As we have noted previously the UV divergences occur in the transverse-momentum integrals 
of partons interacting with a bare cut vertex. To extract this dependence properly it is 
sufficient to separate the perturbative loop from correlation function in question. To this 
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liup 



Figure 6: Modified rules for the discontinuities of the Feynman diagrams. 



end, the latter can be represented in the form of momentum integral in which the integrations 
over the fractional energies of the particles attached to the vertex are removed 



AO 
M(0 



^-X-AC-z) ( i 1 )F(k), 

(2tt) 4 VS ; \ \mi+ ) y h 



where 



F(k) = I d 4 xe ikx (0\^{x)\H,X){H,X\^{0)\0). 



183) 



184) 



In the same way we can easily write down the corresponding expressions for the three- 
particle correlation functions. 

Let us consider, for defmiteness, the fracture function I. Simple calculation of the one- 
loop diagram depicted in Fig. [F] (a) for the 2^2 transition gives in the LLA 



AO 



A 2 



(2tt) 4 v ' J (2tt) 3 
x{-d^(k" - k) lfi (f + m)I{f + m) lv ] 



a 



2?r 2 J (2vr) 4 



^ m I dz w-QH*-'-*) u k „t da , s /„ + *>r- 



2(z" - z) 
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f 



(«) 




Figure 7: One-loop diagrams contributing to the transition kernels of the two-particle 
correlation functions to the two- and three-parton ones. 



1 



[2a" z" + kf} 2 
dz 



[2a" z" + k 



"2l 



2m7 H 



• „ [2aV; + k'j 2 } 
« H s — 



I(z)-M(z) 1 + 2- 



■K-z).' ' (185) 

As long as the logarithmic contribution appears when <C 1 and a /a" <C 1 we 

expand the integrand in powers of these ratios keeping the terms that do produce the 
logarithmic divergence. Similarly, one can evaluate the transition amplitudes of I to the 
three-particle correlation functions Zj (for the diagrammatical representation, see Figs. 
(b,c)): 

a r r 2 I 

J(() A2 = il In A 2 / dzdz'6{(-z)Z 1 {z',z) 
2tt J 

a 



+ 



AO 



A 2 



2vr 



In A / dzdz'6{(- z)Z 2 {z,z') 



(C-z) (z-^') : 
(C-y)(z-^)- 



;i86) 
;i87) 




Figure 8: Contact-type contribution to the evolution equation of the fragmentation function 
T(z). 



Due to the non-quasi-partonic [26] form of the vertex / there exists an additional con- 
tribution to the evolution equation coming from the contact term (Fig. |8|) that results from 
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the cancellation of the propagator adjacent to the quark-gluon and bare cut vertices. As a 
consequence the vertex acquires the three-particle piece 



AO 



A 2 



d A k d A k' 
(2tt) 4 (2tt) 



a 



-Z lp {k', k)igT p {k - k', k)iG{k)6{( -z) + (c.c) 



In A 2 J dz'Z.iz', C) J dz" SKS " J Q m(z", z" - C, z" - C + z>). 



As can be seen Eqs. (|185|) and ( |186| ) possess the IR divergences at z = (. They disappear 
after we account for the virtual radiative corrections (renormalization of the field operators) 



discussed in section 2.6. The net result looks like 



(189) 



Assembling all these contributions we come to the evolution equation for X given below 
by Eq. pH). 

3.4.2 Evolution equations. 




J 




I If 



1 




I 7 





I 1 




I 1 




Figure 9: One- loop (abelian) corrections to the three-point correlation function Z\. 

Now following the procedure just described it is not difficult to construct the closed set 
of the evolution equations by calculating the one-loop diagrams shown on Figs. [|, |T0 



M(0 



a: 
2^ 



dz 



0(C - z)P J 



MM 



M(z) 



(190) 
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Figure 10: The same as in Fig. but for the three-point correlation function Z 2 - 



±(0 = TJ d i^- z A p A^ x{z)+p A^) M{z) 
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+ *(0 
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-^(Z, c - CO + 0(C - 07v7^SvrvP^7T^(s' 
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i (C - 



(193) 



z 
dz' 



C 2 



C 2 z 



+ e° n (C, C - ^)777^a [^(C, *0 - ^(C, CO] 



(C - «0 



±[Z 2 ((,z')-Z 2 ({,(')} 



+ 2^(C,C0 



(C - *0 

where the dot denotes the derivative with respect to the UV cutoff ' = A 2 d/dA 2 and splitting 
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functions are given by the following equations 



Pmm\ z ) 

Pjm{z) 
Pj Zl (z,y) 

Pz±zAz,y) 



z(l-z) 



+ - + 1, 

z 



-1 + -6(z-l) 
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z(l- 


*). 


2 




z(l- 


*). 


2 




z(l- 


*)] 



+ - + 1, 

z 

+ - + T ^—-5(z-l)-ln(l-y), 
z 1 - yz y 



+ - + —^— + 5{z-l) 
z 1 — yz 



- ln(l - y) 



(194) 

(195) 
(196) 

(197) 

(198) 



Now it is an easy task to verify the fulfillment of the equation of motion ( 175 ) for the 
correlation functions as a consistency check of our calculations. By exploiting this relation we 
exclude X from the above set of functions and reduce the system to the basis of independent 
gauge invariant quantities {M and 1Z j}. 

An important note is in order now. As distinguished from the evolution of the structure 
functions, the above Eq. ( |192| ) has the logarithmic dependence on the ratio of the parton 
momentum fractions^. The consequence of its presence is obvious. Taking into account the 
restrictions imposed by Eq. ( |180| ) we can define the moments of the correlation functions in 
the following way 



M n 



We find for two-particle cut vertex 

M n 



a 
2^ 



(S n + S n . 2 ) M n - 



(199) 
(200) 

(201) 



Comparing it with Eqs. ( |115| ) and (|117|) we notice the universality of the evolution kernels 
for the space- and time- like two-particles quasi-partonic cut vertices, i. e. the Gribov-Lipatov 
reciprocity ( |165| ) relation, which looks like 



<y TL -<y Sh 

ln+2 — In 



(202) 



in terms of the corresponding anomalous dimensions, is fulfilled. 

On the other hand, it is impossible to write down the finite system of equations for 
any given physical moment of the three-parton correlation functions as the logarithm of 



9 This is consequence of the fact that while the support properties of the multi-parton recombination 
functions are different from the distributions, i.e. the region of attained momentum fractions is nonsymmetric 
in the former case in contrast to the latter, the perturbative one-loop renormalization of the field operators 
is not sensitive to this. Actually, it is given by the same equations as in the previous discussion of the deep 
inelastic scattering 
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the ratio of the parton momentum fractions in the evolution kernels leads to the infinite 
series of the moments as distinguished from the deep inelastic scattering where the rank of 
the anomalous dimension matrix was finite and increases with the number of the moment. 
This is consequence of essential nonlocality of the cut vertices in the coordinate space since 
even if we start from the local cut vertex it will smeared along the light cone upon the 
renormalization. Therefore, it is not possible to solve the system of equations successively in 
terms of moments as well as we do not succeed in solving it analytically in a general form. 
However, in the next section when dealing with the QCD evolution we will find that the 
system of coupled equations ( |192| ) and (|193|) can be reduced to the single equation in the 
multicolour limit and its solution can be found analytically. 



3.5 Non-abelian evolution. 

In the QCD case we should add the diagrams with triple-boson interaction vertex (see Figs. 



IT, O) and gluon self-energy insertions (not shown). 




Figure 11: Non-abelian radiative corrections to the evolution kernels of the fragmentation 
function Z\. 




Figure 12: The same as in Fig. [11] but for the fragmentation function Z^. 
Gathering these contributions together with equations obtained in the previous section 
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(with colour group factors accounted for properly), we come to the final result 



a; 



2tt n 

+ / ~ z) 



©Si (co c - OtMO - °(C)t77^m(( - a 



C(C-C')' 



+ f (- 2 4/o 1a,Zi(c '-" (c - ;) ' 2)+ ( 

+ y ^ 



(* - C + CO cc 



-{ + {', z) 



-c F ( e^c, C - O^^W, + 0(0 C - o) 



+ Cf-^ e° u (C,C- CC'-C + ^O 



^ (C - C + *0 



Ca ( i /-I H a H i\ z \ z + C ~~ C) ~ / / >n 



r\0 i H a' /-\ (C C )(C + z ) ~ i i /-\ i ,^0 ( /-i /-i i\ (C + C) 

+ ©ii(C,C - O777- ' " 



0(C) 



(<;>-z')(C-z') 

i (C-C' + z' 



+ ©i\(C, C'-z'y^f^iz'x) 



.(z'-(') 



(z> 



z^z'-CX-C) 



-2e° u (C,C-z') [z, (z',o-z 1 (C, ()} 



(203) 



MC,0 = ^\-CfO(c-0- i m(0 



c 



+ 



dC 1 



'c c ,s 



C ^* ( *i^U(*-C,*) 



+ ^ (-2^ jf 1 d^(«, c - «(c - z)) + 



c 2 



C F 6((-C) C -^Z 1 (z>,C) 



+ c 



c. 



o(C - CO 



, N (C-C0(C-C' + *0 



CV 



>-c+co c 



Zi(z',C) 



z -^)z 2 (z,z-( + a 



+ e^c, C - z') j^^iMC, *0 - 3»(C, CO] 



+ ©i\(C - C, C - ^)|^[^(C, *0 - ^(C, CO] 



* (c - CO ( i 



2 c,) (c-^o l(C'-*0 1 c 



—r + - Z^Z'X') 



(204) 



41 



where 
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These equations should be supplemented by the equation for the mass dependent correlator 
Ai which differs from its abelian analogue (|190|) only by the Casimir operator Cp- 



3.6 Asymptotic solution of the evolution equations. 

One can easily observe the significant reduction of the above evolution equations if we neglect 
the terms in the kernel of the order of magnitude 0(1/N%). In this case an additional three- 
parton correlator Z 2 ~ (OliptplH, X)(H, X|5 J -|0), which appears only through the radiative 
corrections, decouples from the evolution equation for Z\. Therefore, discarding the quark 
mass cut vertex we obtain homogeneous equation which governs the (^-dependence of the 
three-parton correlation function Z\. 

The situation has the closer similarity with phenomenon found in the evolution equations 
for chiral-even and -odd distribution functions discussed in the first part of this paper [^, 0, 
38[ where in the multicolour limit (iV c — ► 00) there was a very important simplification as the 
evolution kernels have been vanishing for contributions with interchanged order of partons 
on the light cone, i.e. the momentum fraction carried by gluon in the matrix element of 
quark-gluon correlator varies only among the quark ones and does not exceed the latters. 
This property allowed to find the solution of simplified equations exactly in the nonlocal 
form. In the present case the decoupling of Z\ has the same consequences. 

In the large- N c limit the RG equation takes the form 



An J z 



z)1C{z,z!^,C)Zi{z\z) 



and the evolution kernel is given by the following expression: 
1 



N, 



■1C(z,z',(,C) 




(206) 



(207) 
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Inspired by our knowledge acquired from the previous study of the twist-3 structure 
functions we are able to check that Eq. ( |176| ) (with M. = 0) satisfies the ladder-type 
evolution equation with the following splitting function: 



§(>-§) 



Thus, for the moments we obtain the following solution of the RG equation (Q > Qq) 
and the corresponding anomalous dimensions equal 



(209) 



NC °7„ = N c [-2^{n - 1) - 2 lE - + 1} , (210) 

as usual /3 — f ^/ — y^A- 

As we have previously mentioned, there exists an equation which states that in the 
leading log approximation the time-like (TL) and space-like (SL) kernels corresponding to 
the twist-2 parton densities are directly related by the Gribov-Lipatov equations (|165|) or 
( |202| ). Comparing the result given by Eq. ( |210| ) with the large- iV c anomalous dimensions 



( 163|) we see the absence of the universality of the corresponding twist-3 evolution kernels, 



i.e. the Gribov-Lipatov reciprocity is violated. 

3.7 Generalization to other fragmentation functions. 

Now we can proceed further and demonstrate that the evolution kernels for the time-like 
two-quark densities can directly be found from their space-like analogues by exploiting the 
particular form of the evolution equations given by Eqs. ( |161| ) and (|162|) . Since the analytic 



structure of the uncut diagram (see Fig. |13|) is completely characterized by the integral 
representation of the 9-function given by Eq. (|69"D , we can just take its particular disconti- 
nuities, using the usual Cutkosky rules supplied with appropriate theta-function specifying 
the positivity of the energy flow from the right- to the left-hand side of the cut, in order to 
obtain the corresponding time-like kernel. Since the observed particle is always in the final 
state for the fragmentation process, we are restricted to the single cut0 across the horizon- 
tal rank of the ladder diagram in Fig. Namely, using the integral representation of the 
corresponding step function 0° l5 we have 

Q0 (xx _ g)= f 00 da 1 6{x - g) 

111 ' P) J-a a 2m[ax-l + i0][a(x-P)-l + i0] (3 ' 1 ' 

The self-energy insertions are not affected by the cut since it does not cross the corresponding 
lines. Taking into account different kinematic definitions of the correlation functions in the 
space- and time-like regions^ H , we are able to find the kernels. It is easy to verify that the 
evolution kernels constructed for the time-like twist-2 cut vertices using this recipe coincide 
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Figure 13: One-loop ladder-type diagram for the two-particle evolution kernel. 



with the known results. In the same way, we may obtain the above equation ([208 ) from Eq. 

Since there exist fragmentation functions corresponding to each distribution, apart from 
the specific ones appearing from the final state interaction, we are in a position to find 
large-iV c anomalous dimensions, which govern their Q 2 -dependence, from the results (|162|) 
and Qiuip . Namely, the genuine twist-3 contributions Hf~ 3 (PCO) and G^ 3 (PCE) to the 
corresponding fragmentation functions 

Hl(() — 7 / ^e iX <(0\ia + ^(\n)\H,X)(H,X\i>(0)\U), (212) 

MO = \j ^e iX <(0\ 1±1 ^(Xn)\H,X){H,X\m\0)- (213) 

after subtracting out the twist-2 piece |28| obey the evolution equation ( [209| ) with the fol- 
lowing anomalous dimensions: 

PC °7„ = iV c {-2V(n-l)-2 T£ + -l I + l}, (214) 

PCE 7„ = iV c {-2^n-l)-2 7E --^ + ~}. (215) 

Of course, it is a trivial task to invert the moments and to find the DGLAP kernels them- 
selves. 

To summarize this section, we have found that in the multicolour limit of QCD the twist- 
3 fragmentation functions obey the ladder-type evolution equations and the corresponding 
anomalous dimensions are known analytically. The Gribov-Lipatov reciprocity is not the 
property of the twist-3 distributions but it is strongly violated already in the LLA of the 
perturbation theory. 



4 Discussion and conclusion. 

We review above the approach to an analysis of the logarithmic violation of the Bjorken 
scaling in the twist-3 distribution and fragmentation functions of the nucleon. It consists in 

10 For a given graph the possible cuts correspond to the possible final states. 
n See, for instance Eqs. @ and ( |l70| ). 
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the studying of the one-loop renormalization of the multi-parton correlators which explicitly 
involve the gluon degrees of freedom and further reconstruction of the evolution equation for 
them in the LLA using the renormalization group invariance. For these purposes, we have 
used the techniques, which employ the light-like gauge for the gluon field. The physically 
transparent picture which appears in this gauge (which is an essential ingredient of our 
method) makes the calculations simple. Accepting different prescriptions on the spurious IR 
pole in the gluon propagator, we were able to verify that they do lead to the same results. 
We present an exact leading-order evolution for the correlators in the light-cone fraction as 
well as in the light-cone position representations, which display the complementary aspects 
of the factorization, and establish the bridge between different formulations of the QCD 
evolution. 

From the calculational point of view the momentum space technique is much easier to 
treat. However, the coordinate space makes the involved symmetries apparent and, as a 
by-product, diagonalization of evolution kernels is easy to handle. The complicated form of 
exact master equations for the twist-3 functions compels one to look for the approximation 
which could work with reasonable accuracy. The solution of this problem has been found in 
the fact that in the multicolour limit of QCD they are reduced to the ladder-type evolution 
equations which generally have very good precision, at the level of few per cent. 

Comparing the analytical expressions for the anomalous dimensions of the twist-3 struc- 
ture and fragmentation functions we observe that the Gribov-Lipatov reciprocity relation, 
valid for low-twist parton densities in the LLA, is broken already in the leading-order of 
perturbation theory. 

The result we have discussed here are important from the theoretical point of view since 
they enrich the theory of the higher-twist effects in the hadron reactions as well as for 
phenomenology and could be used to analyze the experimental data when these become 
available. 
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